Toda Hierarchy with Indefinite Metric 

Yuji KODAMA* and Jian YE+ 

Department of Mathematics, Ohio State University, 
Columbus, OH 43210, USA 

May 15, 1995 



Abstract 

We consider a generalization of the full symmetric Toda hierarchy 
where the matrix L of the Lax pair is given by L = LS, with a full 
symmetric matrix L and a nondegenerate diagonal matrix S. The key 
feature of the hierarchy is that the inverse scattering data includes a 
class of noncompact groups of matrices, such as 0(p,q). We give an 
explicit formula for the solution to the initial value problem of this 
hierarchy. The formula is obtained by generalizing the orthogonaliza- 
tion procedure of Szego, or the QR factorization method of Symes. 
The behaviors of the solutions are also studied. Generically, there are 
two types of solutions, having either sorting property or blowing up 
to infinity in finite time. The r-function structure for the tridiagonal 
hierarchy is also studied. 
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1 Introduction 



The finite non-periodic Toda lattice hierarchy can be written in the Lax form 
Hi) Hi with variables t := (£1,^2, • • •)> 



n 



1,2, 



(1-1) 



where L is an N x N symmetric "tridiagonal" matrix with real entries, 

/ 01 61 ... \ 
61 a 2 b 2 ... 

: '•■ : 

cln-i &jv-i 

fejv-i oat y 

and _B n is the skew symmetric matrix defined by 
B n = Y[L n := (L n ) >Q -(L^ 



V 



;i.2) 



'<o 



;i.3) 



Here (L n ) >0 , <0 x denotes the strictly upper (lower) triangular part of L n . (We 
formally write an infinite number of flows in ( |1.1|) , even though there are at 
most N independent flows.) In particular, (|1 . 1| ) for n=l with t := t\ has the 
form 

da k 



dt 



-— — Ok\a k+ i — at ■ 
at 



;i-4) 



;i.s) 



This system describes a hamiltonian system of N particles on a line inter- 
acting pairwise with exponential forces. The hamiltonian for the system is 
given by 



A? 



N-l 



H = 9 J2 vl + £ exp(x fc - Xk+i] 
/ fe=i fc=i 



where the canonical variables (xk, yu) are related to (a^, 6^) by 



3/k 
' 2 : 



;i.6) 



;i-7) 



and 



1 (x k 
- cxp 



Xk+l 



;i.8) 



The r-functions were introduced in || to study the Toda equation (|1.4j) 
and ( |1.5|) . Writing a^s with the r-function Tj's and t = ti, 



hi can be expressed as 



19, ^ 

2 at Tj_i 



2 lr i+1 Ti_i 



4 r 2 



;i.9) 



(1.10) 



Then (L.4) and (|1.5|) become 



1 <9 2 r i+ iTj_i 
4 at 1 rf 



[1.11) 



These r-functions Tj have a simple structure, that is, a symmetric wroskian 
given by 0, 



fll #2 

01 #2 #3 

ft-i ft • • • 



ft-i 
9i 



(1.12) 



where the entries g n are defined by g n := (l/2)dg/dt n , with a function g 
satisfying the linear equations 



1 dg 1 <9™# 
2d£ ~ 2™~dF' 



(1.13) 



These r-fuctions are also shown to be positive definite in ||7|. It is a key prop- 
erty of their relation to the moment problem of Hamburger. The equation 
( |1 . 1 1| ) is also expressed by the well-known Hirota bilinear form 



Dt T i ' T i — $>Ti + lTi-l 



;i.i4) 
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where the Hirota derivative is defined by 



(D t f.g)(t) := ±f(t + s)g(t-s)\ a=0 

~ <™> 

There have been extensive studies on the hierarchy ( |1.1Q and its general- 
izations. One of them is to extend L from "tridiagonal" to "full symmetric" . 
This more general system, which we call "full symmetric Toda hierarchy", 
was shown by Deift et. al. in JE] to be a completely integrable hamiltonian 
system. The inverse scattering scheme for ( |1.1[ ) with L being any symmetric 
matrix consists of two linear equations, 

L$ = <M , (1.16) 

= *.*. (1-17) 

where $ is the orthogonal eigenmatrix of L, and A is diag (Ai, ■ • ■ , A at)- The 
hierarchy (|1.1|) then results as the compatibility of these equations with 
dA/dt n = (iso-spectral deformation). In [|J, Kodama and McLaughlin 
solved the initial value problem of ( 1.16Q and ( 1.17 ) using the "orthonor- 



malization method", and derived an explicit formula of the solution in a 
determinant form. They also showed that the generic solution assumes 
the "sorting property". Here the sorting property means that L(t\) — > 
A = diag (\\, ■ ■ ■ , \n) as t\ — > oo, with the eigenvalues being ordered by 
Ai > A2 > ■ • • > Ajv- 

The hierarchy we consider in this paper is for L = LS, where L is full 
symmetric and S is a constant diagonal matrix, S = diag (s\, ■ ■ ■ , sn) with 
nonzero "real" entries. All the entries of L are assumed to be real, unless 
otherwise stated. The hierarchy is also defined as the Lax form (|1.1|) with 
( |1.3| ). In particular, for the case of tridiagonal L with n=l, we have 

^=2(s k bt-s k ^bl_ 1 ), (1.18) 
^ = b k (s k+1 a k +i - s k a k ). (1-19) 
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In analogue to the nonperiodic Toda equation with ( |1.6| ), the hamiltonain H 
for the above equations is given by 

i N N-l 

H = 9 Vk + SkSk+i exp(x k - Xk+i), (1-20) 
z k=i fc=i 



with the change of variables 



SkCik = -y, (1-21) 



and 



SkSk+ibk = ^ ex P I ^ ) • \ 1 - 22 ) 



Note from ( |1.20| ) that H includes some attractive forces, thereby is not pos- 



itive definite. One then expects a " blowing up" in solutions. 

We study the initial value problem of the hierarchy by the "inverse scat- 
tering method". The inverse scattering scheme is also given by (|1.16|) and 
( |1.17 ), where the eigenmatrix $ is now normalized to satisfy 



QS- 1 ^ =S-\ $ T S$ =S. (1.23) 

In the case of the full symmetric Toda hierarchy, where S is the identity 
matrix, $ is given by an orthogonal matrix, $ e 0(N). As a special case of 



1.23|) , we have $ G 0(p, q) for S = diag(l, —1, • • • , —1), with p + q = 

N. From the first quation in ( |1.23| ), we define an inner product where S^ 1 
gives an indefinite metric for S being not positive definite. For this reason, we 
call the hierarchy considered in this paper the "Toda hierarchy with indefinite 
metric" . 

The content of this paper is as follows: We start with a preliminary in 
Section 2 to give some background information on the Toda hierarchy and 
the inverse scattering scheme. First we show the compatibility of the flows 
in ( |1 . 1|) . Then we show that the eigenmatrix of L can be chosen to satisfy 
( |1 . 23|) , and also (|1.23|) is invariant under the flows generated by ( |1.17|) . 

In Section 3, we show that the hierarchy can be solved using the orthonor- 
malization method with respect to the indefinite metric. The solution of the 
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hierarchy turns out to be given by the same solution formula in PJ with a 
little modification. 



In Section 4, we discuss the behaviors of the solutions. Due to the in- 
definiteness of the metric defined by S 1-1 , the solution has richer behaviors 
than the full symmetric case. Generically, in addition to the sorting property, 
there are solutions blowing up to infinity in finite time. The behaviors of the 
solutions are characterized by S and the initial conditions, i.e., the eigen- 
values Aj, i = 1, • • • , N of L and the initial eigenmatrix $°. In this paper, 
we obtain the following main results: If S is positive definite, then generic 
solutions have the sorting property (Theorem 2). If some eigenvalues of L 
are not real, or $° is not real, then generic solutions blow up to infinity in 
finite time (Theorem 3). 

In Section 5, we illustrate these results with explicit examples. 

In Section 6, we study r-function structures of the L hierarchy. In the 
case of tridiagonal L, one can also introduce r-functions in the same form 
as the symmetric wroskians (|1 . 1 2|) . However, our r-functions are no longer 
positive definite as a result of the indefinite metric in our hierarchy. We also 
find a "bilinear identity" generating relations among Tj's of ( |1.12j ). From the 
bilinear identity, we then derive a hierarchy written in Hirota's bilinear form 
for the r-functions including (|1.14j ) as its first member. 



2 Preliminary 

The hierarchy considered in this paper is defined as 

d ~ 



dt n L 



B n , L 



n = l,2,..- (2.1) 



where L = LS with L, a full symmetric matrix, and S, a nondegenerate 
diagonal matrix, i.e., a symmetrizable matrix in the sense of Kac ||, and 
B n := Y\ a L n = (L n ) >0 — (L") <0 . One should note that among the infinite 
number of flows in the hierarchy, only the first N flows are independent, and 
the higher order flows are related to these N flows through P(L) = 0, the 
characteristic polynomial of L. 
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Let us first show the compatibility of flows in ( |2.1| ), that is, d 2 L/dt m dt n = 
d 2 L/dt n dt m . This can be shown from the "zero curvature" condition, i.e. 

J^T ~ = ^ m ' ^n] ■ ( 2 - 2 ) 



Proposition 1 The flows in ( \2.1\) commute with each other. 

Proof. It suffices to show the zero curvature condition fl2.2|) . From Q2.1|) , 
we have 

df n df m 

— - — = lB m ,L n }-lB n ,L m }, (2-3) 

(Jim (JTi n 

which can be also written as 

[B m , L n ] - [B n , L m ] = [B n , B m ] + [B m + L m , B n + L n \. (2.4) 

Note that B m +L m is an upper triangular matrix. Taking the lower triangular 
projection of (|2.3|) , we have 



d(L n ) <0 d(L m ) <0 
dt dt 



Similarly, 



d(L n ) >0 d(L 



>0 



([B n ,B m }) <0 . (2.5) 



-([B n ,B m )) >0 . (2.6) 



dt dt 
Subtracting ( |2.6| ) from ( |2.5| ), we get 

o7~" = {[Bm ; -Bn])>0 + ([-^m ) B n ]) <0 . (2.7) 

Since [JB m , £? n ] is skew-symmetric, diag([B m ,B n ]) = 0. This completes the 
proof. □ 

Remark 1 Since the specific form of L is not used in this proof, Propo- 
sition 1 is valid for an arbitrary matrix. 

In the case of L being full symmetric, L can be diagonalized by an orthog- 
onal matrix, and the flows defined by (|1 . 1|) are compatible with the choice 
of skew-symmetric B n by ( |1.3|) . To set up the inverse scattering scheme for 
our system, we need the following lemma from linear algebra: 
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Lemma 1 Let L be a matrix given by L = LS , where L is symmetric, 
S = diag (s\, ■ ■ ■ , sn) nondegerate. Suppose that all the eigenvalues of L 
be distinct. Then L can be diagonalized by a matrix $ satisying i.e., 

qs- 1 ^ = s- 1 , $ T S$ = S. 

Proof. Let (■,■) be the usual euclidean inner product, i.e. (x, y) : = 
J2iLiXiVi, for x = (xi,---,x N ) and y = (yi,---,y N ). Let 0(A;) be the 
eigenvector of L corresponding to the eigenvalue Aj, i.e., L(f>(\i) = Aj0(Aj). 
Calculating (S(f>(\j), L0(Ai)), we have 

(50(A,),Z0(A J )) = A 4 (50(A i ),0(A 4 )) 

= (LS'0(A i ),50(A i )) = A,-(50(A i ),0(A i )), 

which leads to 

(A J -A,)(50(A i ),0(A l )) = O. 

By the assumption, Aj 7^ Xj for % 7^ j, we obtain (S<f)(\j), 0(Aj)) = 0, that is, 
$ T 5$ is a diagonal matrix. It is also nondegenerate. We then normalize $ 
to satisfy $ T S$ = S. We multiply S' _1< I )T on the both sides from the right, 
get ^S^S- 1 ^ = $ T . This leads to ^S' 1 ^ = S' 1 . □ 

Remark 2 With the normalization, the eigenmatrix $ becomes complex 
in general, even in the case that all the eigenvalues are real. This is simply 
due to a case where the sign of Y,k=i s k4>l{\) m ® T S<& differs from that of 

Si. 



We now show that the choice of B n in (|1.17|) is compatible with (|1.23|) , 
that is: 

Lemma 2 Eqs (\1 . 23j ) are invariant under the flows generated by B n = 

I\ a L n . 

Proof. First, we have 

dt n dt n dt n 

= <S> T (BTS + SB n )<5> . 

So all we need to show is B^S + SB n = 0. Since S is diagonal, it commutes 
with the anti-symmetric projection ]J a} that is, 

B n :=(L n ) >0 -(L n ) <0 = (L^LSU-iL^LS)^ 

= (L^LUS-iL^L)^. 
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This leads to 

BlS + SB n = S[{L n - l L) >0 - {L n - X L) <Q \ T S 
+S[(L n - 1 L) >0 -(L n - 1 L) <0 }S. 

Since L n ~ 1 L is symmetric, (L n ~ 1 L) >0 — (I/ n ~ 1 L) <0 is skew-symmetric. Thus, 
we have B*S + SB n = 0. □ 

The eigenmatrix $ consists of the eigenvectors of L, say (f)(Xk) = (0i(Afc), 
•..,^(A fc )) T forA; = l,2,...,iV, 

$ = [0(Ai), ••• , <K\ N )] = MWi^n • (2-8) 
Then ( 1.23|) give the "orthogonality" relations 

JV 

E*fcVi(A*)0i(A Jfe ) = < yysr 1 > (2-9) 
fe=i 

^ s fc 0fc(A i )0fc(A j ) = 5ySj . (2.10) 

k=l 



With ( p.9[) , we now define an inner product < •, • > for two function / and 
g of A as 



N 



</, 2 7(A*)s(A fc ), (2.11) 

fc=i 

which we write as < fg > in the sequel. The metric in the inner product is 
given by 

N 

dor(A) = $>^(A- \ k )d\ , (2.12) 
fe=i 

which leads to an indefinite metric due to a choice of negative entries in 
S. The entries of L n are then expressed by 

~af := (Z n ) .. = Sj < X n Mj > • (2.13) 
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3 Inverse scattering method 



In this section, we solve the initial value problem of the hierarchy ( |2.1|) by 
using the inverse scattering method. Namely, we first solve the time evo- 
lution for the eigenmatrix $(t) of L with the initial matrix $(0) := $° = 
[4 ) i(^j)}i<i,j<N obtained from the eigenvalue problem L(0)$° = $°A, and 
then find the solution L(t) through ( [2.13| ) with $(t). Here we call the eigen- 



values A*., k 



N and the eigenmatrix $ the "scattering data" . 



Let us first note B n = L n — diag(L n ) — 2(L n ) <0 , and using L<3> = $A, we 
write ( |1.17[ ) as 

° \diag(L n ) +2(L n ) <0 



dt n 



Then the equations for the first row vector 0i(Afc,t), k — 1, 



(3.1) 
N, in $(t) 



arc 



= (AS - -i < A»tf(A) >)0!(A fc ), 



which can be readily solved in the form 0i(Afc, t) 

^ 1 (A fc ,t)=0?(A fc )e« A - t ), 



(3.2) 



x / S i<V?(A,t)> 



with 



(3.3) 

where e^ Afe,t ' ) := Z^^Li A£t n . The equations for the second row vector ^(Afc, t), 
k — 1, ■ ■ ■ , N, are 



= (AS - s 2 < A n ^(A) >) </> 2 (A fc ) 
-2s x < A> 2 (A)0 1 (A) > 0x(A fc ) . 



The solution to 
with 



can be also found in the form 2 (Afc,t) 



(3.4) 

■02(Afc,t) 
y/ S2<1pl{\,t)> 



^ 2 (A fc ,t) = ^(A*)^) - Sl < 0°(A)e^ A ' t Vi(A,t) > &(A*,t). (3.5) 
In general, the ith row vector of $ in ( |3.1| ) satisfies 

d(t>i{\ k ) _ f xn „ ^ xn ,2/ 



A£- Si <A"0, 2 (A)> WA fc ) 



i-l 



-25><A n &(A)0,(A)>&(A fe ). 



(3.6) 
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Then the solution of ( |3.6| ) is expected to have the form (fii(\k, t) — ., , ^ 

with 

Uh, t) = $(\ k )e^ - £ Sl < 0O(A)e« A ' t VKA, t) > <MA fc , t) . (3.7) 

As we will see below, this is indeed the case. The above procedure sug- 
gests that the eigenmatrix $ can be constructed by the "orthonormalization 
procedure" on the row vectors [$ (Xk)^ Xk '^]i<k<N , i = 1, ■ ■ ■ , N. 

Let us now give a more systematic presentation of the above procedure, 
which is the essence of the approach in || . Because of the simple structure 
of ipi(\k) in ( |3.7|) , we first write 



$ = Tty , (3.8) 

where * := [i>i(\j)]i<i,j<N and 

T = drag \{s x < i>\ >)-^ 2 , • ■ ■ , (s N < ^ >)" 1/2 " . 



Note that ( |3.8j ) can be regarded as a gauge transformation and includes a 
freedom in the choice of ip. Namely, (|3.8j ) is invariant under the scaling of ipi, 
i.e., fa -> /i(t)V»i, with {/i}^! arbitrary functions oft. With M, ( PTTEQ 
and ( |1.17|) become 

(t _1 LT) * = A , (3.9) 

A* = (T-'^.-fAtogr)*. (3.10) 

Writing as in ( |3.1| ), i.e. 

T- l B n T) = -2 (r- l L n T) <Q + (T^U'T) - diag (L r 



( |3.10p gives 

dtp 
dt n 



-2 (T- 1 L n T) <Q i/> + A> - ^ (L") + A logT ^) ^ . ( 3.ii) 
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We here observe that ( |3.11| ) can be split into two sets of equations by fixing 
the guage freedom in the determination of ip. In the components, these are 



dipt 
dt n 



Id o ~(n) 

2dt n 6 Yl 



(3.12) 



(3.13) 



It is easy to check that (|3.12|) implies ( |3.13|) . It is also immediate from ( |3.12j ) 
that we have: 

Proposition 2 The solution of ( \3.13{ ) can be written in the form of separa- 
tion of variables, 



V>(A,t) = A(t)<f(\)e 



(3.14) 



where A(t) is a lower triangular matrix with diag[A(t)] = A(t = 0) = 1^, 
the N x N identity matrix, and 0°(A) = 0(A, 0). 

Note that we have chosen the initial conditions of ip(\,t) to coincide with 
0(A,t), i.e. 4>(\,0) = 0°(A), and thus s { < ipiipj > (t = 0) = s { < >= 
5ij. As a direct consequence of this proposition, and the orthogonality of the 
eigenvectors, ( ^9|) , i.e. < ipiipj >= for i ^ j, we have: 



Corollary 1 (Orthogonality): For each i G {2, • ■ ■ , N}, we have for all t 
with t m G R, 



(3.15) 



fc=i 



for j = 1,2, 



2-1. 



Now we obtain the formula for the eigenvectors of L in terms of the initial 
data {0°(A)}i<j<Ar: 

Theorem 1 The solutions ipi(\,t) of ( \3.12[ ) are given by 

SlCn . . . S\Cu 



<MA,t) 



e «(A,t) 
A-i(t) 



0?(A) 



#(A) 



(3.16) 
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where %(t) =< 0°0°e 2 ^ A ' t - ) >, and A(t) is the determinant of the k x k 
matrix with entries SjCy(t) ; i.e., 



D k {t) 



(siCij (t 



l<i,j<k 

(Note here that SjCjj(O) = 5^ and A(0) = lj 

Proof. From equation ( |3.15| ) with ( |3.14j) , we have 

s t £ Aih(t) < <P k<py^ X,t) >= , for 1 < i < i - 1 . 
fc=i 



(3.17) 



(3.18) 



Solving (|3.18|) for Aik with Ajj = 1, we obtain 

AM : 



A-i(t) 



(3.19) 



where A) 00 := 1, and Df_ 1 (t) is the determinant A-iOO in the form 
( |3.17|) with the replacement of the kth column (sxCifc, . . . , Sj_xQ-i /c) T by 
(sicij, . . . , j) T . From fl3.14| ), we then have 



o 

ik<Pk 



fc=l 



e £(A,t) 
A^7 



A- j A-i 



X -0? 



SlCii_i 



+0° 



sicn • • • s x cn_i 



+ 



which is just (|3.16| ). □ 

We note from ( |3.16| ) that < ipf > can be expressed by A, 

D, 



< Vi > - 



SiA-i 



(3.20) 
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This yields the formulae for the normalized eigenfunctions 



i(A,t) 



MM) 



A(t)A-i(t) 



S1C11 



Si-iCi-n 
0?(A) 



Si— lCj— \i 

$(A) 



(3.21) 



With the formula ( |3.21p , we now have the solution fl2.13| ) of the inverse 
scattering problem ( |1.16| ) and (|1 . 1 7|) . 



The above derivation of fl3.21| ) is the same as the orthogonalization pro- 
cedure of Szego 0, which is equivalent to the Gram - Schmidt orthogonal- 
ization as observed in [Q, except here the orthogonalization is with respect 
to an indefinite metric. Indeed, from the form of </>,■ = , ^ „ . we see that 



( |3.14|) expresses the relations 



4>i(\,t) G span |0°(A)< 



,€(A,t) 



(3.22) 



where we are viewing {0j(A, t)}^ as a collection of "orthogonal" functions 
defined on the eigenvalues of the matrix L. The relations ( |3.22|) together 

s" 1 ^-, imply that {<fo(A,t)} i=1 are 



< 



with the orthogonality (J27 

obtained by an orthogonalization of the sequence {0°(A)e^ A,t )}. , and hence 
one may obtain ( |3.21 ) from the classical formulae as in [Q] . If in the orthogo- 
nalization procedure for 4>i(X), we replace ^4>i(\)e^ x,t \ • ■ • , 0°_ 1 (A)e^^ A,t )| by 
{(pi (A, t), • ■ • , i _ 1 (A, t)}, then we get the solutions given by (|3. 



Remark 3. The above method is a natural generalization of the QR 
factorization method of Symes || . In the case of L being symmetric, the QR 
factorization method is given as follows: First factorize e L(0) * = Q(t)R(t), 
where Q(t) G SO(N) with Q(0) = In, and R(t) is a lower triangular matrix. 
Then L(t) is given by L{t) = <5 -1 (t)L(0)(5(i). In the case of our L, the same 
procedure applies where Q(t) satisfies ( |1.23|) , i.e. $(t) = Q^it)^ . 

Remark 4. The above method applies directly to the following equation, 



Mi 



Bi, L + f(L), 
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where a "pumping" term f(L) is given by an entire function of L. Due to 
the presence of f(L), the eigenvalues of L(t) are no longer time- independent, 
In fact, we have 

^ = /(A), 

and the solution formula fl3.21| ) still holds with the replacement of e xtl by 
exp (^Jq 1 A(t)dt), and U — except i — 1. 

4 Behaviors of the solutions 

Here we consider the first equation of the hierarchy, that is, we have only 
one time ti, which we denote by t. As a consequence of the explicit formula 
(|3.21 ), we can now study the behaviors of the solutions. First we note: 



Proposition 3 D i} i = 1, 2, • • • , N , in are real functions. 



Proof. In the construction of the solutions the "gauge" T is fixed by 
( |3~T3D . In terms of A, (H) is 



^ = 2^ l0g A^" (4 ' 1} 

Note that D = 1, -Dj(O) = 1 and are real functions. Then we see by 
induction that all D; are real functions. □ 



In (|4~T1 ), suppose Di(t ) = for some finite t and some i. Then if L(t ) 
is a finite matrix, Di_i(t ) must be also 0. By induction, D 1 (t ) = 0, but 
D (t) = 1, this forces a n to be infinite, which is a contradiction. So we have: 

Proposition 4 Suppose Dj(io) = for some t < oo and some z ; then L(t) 
blows up to infinity at to. 

We note that Di, i = 1, 2, ■ • ■ , N, are zth principal minors of the product 
of matrices SQeS^&J , where <3> e is 

/ e Al *0i(Ai) • • • e A -*0i(Av) \ 

v e^%(Ai) ••• e A - t ^(A Ar ) / 
Then from the Cauchy-Binet theorem, we have: 
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Proposition 5 Di, i = 1, 2, • • • , N can be expressed as 



,(4-2) 



where Jn represents all possible combinations for 1 < j\ < • • • < j$ < N . In 
particular D (t) = 1, and D N (t) = exp(2X^i \t). 

This proposition is very useful to study the asympototic behavior of Di for 
large t. To study the time evolution of Di(t), we need the information on the 
scattering data, that is, the eigenvalues A& and the eigenmatrix <3>°. Here we 
have from the linear algebra: 

Lemma 3 Let L be a matrix given by L = LS , where L is symmetric, 
S = diag (si, ■ ■ ■ , sn) positive definite. Then all the eigenvalues of L are 
real, and the initial eigenmatrix $° of L satisfying ( \1 . 23j ) is also real. 

Proof. Let <j)(\i) be an eigenvector of L corresponding to Aj, i.e., 

= Xnf>(\). (4.3) 

Multiplying by the adjoint of S<j)(Xi), i.e., (j) T (Xi)S to the left, we get 

T (A,)5L50(A 4 ) = A l T (A,)50(A l ). (4.4) 



On the other hand, if we take the adjoint of ( |4.3j ), and then multiply S(f)(\i) 
to the right, we get 

T (A,)5L50(A 4 ) = A^ T (A l )50(A l ). (4.5) 



Subtracting (|J) from flPp , we have (A; - A;)0 T (Ai)S0(Ai) = 0. Since S is 
positive definite, this forces A^ — Aj = 0, that is, A« is real. 

As for the normalization of 0(Aj), since both T (Aj)5'0(Aj) and Si are 
positive, we only need to multiply 0(A») by some "real" factor. So 0(Aj) 
remains real after normalization. This verifies the assertion of the lemma. □ 

We now obtain: 
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Theorem 2 Let the eigenvalues of L be ordered as Ai > A2 > •■• > Xn- 
Suppose that S is positive definite and <iet$° 7^ for n = 1, . . . , N, where 
$° is the nth principal minor of $°. Then as t — > 00 , the eigenf unctions 
(f)i(\j,t) satisfy 



(f)i(\j,t) — > Sij x sgn (det^fj sgn (detQ®^ 
which implies L(t) —>■ diag (Ai, . . . , Xn) from 



(4.6) 



Proof. Using Lemma 3, we see that all the terms in the sum ( }4.2| ) are pos- 
itive, thereby D^s are positive for all t. From the ordering in the eigenvalues, 
we see that the leading order for Di is given by 

2 



0?(Ai) ... #(A0 



(4.7) 



#(Ai) ••• 0°(A.) 
where we have assumed detQ® 7^ 0. From (ft.21|) and (|4.7|), as t — > 00, 



0n(A;t) 



X 



|de*$° | Ide^il exp [(2 E^ 1 A, + A n ) t 

SiCn(t) •■■ SiCi n (t) 
$11— lCln— 1(0 ' ' ' Sn—lCn—ln 

0?(A) ••• 0°(A) 



x 



(4.8) 



The dominant term in the determinant gives 



0?(A) 



^(A.J 

<(A^ 
0°(A) 



0?(Ai) 



<(Ax) 



x 



0?(A n _O ... 0°(A„_!) 
0?(A) • • • 0°(A) 

p»-l 



(4.9) 
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where a({£j}™ =1 ) is the signature of the permutation P n _i for 1 < lj < n — 1. 



1. 



We note that the determinant in (|4.9| ) is zero for A 
The ordering Ai > • • • > Ajy implies the result. □ 

Theorem 2 shows if S is positive definite, then generic solutions have the 
"sorting property". It is a natural generalization of Theorem 2 in |4j, where 
S = In- Next theorem provides sufficient conditions for the solutions to blow 
up to infinity in finite time. 

Theorem 3 Suppose some eigenvalues of L are not real, or <3>° is not real, 
and detQ° n ^ for n = 1, ■ • • , N. Then L(t) blows up to infinity in finite 
time. 

Proof. We have two cases to consider: 

a). All the eigenvalues are real, but $° is complex. 
From Remark 1, each column 0°(Aj) of $° is either pure imaginary or real. 
In this case, let k be the first column such that <p°(\ k ) is pure imaginary. 
With the ordering Ai 
of D k is still given by 
imaginary, while </>°(A, 



> 



> Aat, the leading order in the expansion ( |4.2j ) 
We also assume det<&° k ^ 0. Since 4>°(X k ) is pure 



k — 1, are real, we have 



0\2 



0?(A fc ) 



< 



Being negative of the leading order in Dk implies lim t ^oo Dk(t) = — oo. Note 
D k (0) = 1, so there is some t < oo, such that D k (t ) = 0. Then Proposition 
4 results that L(t) blows up to infinity in finite time to- 

b). Some eigenvalues are not real. 
We order the eigenvalues of L by their real parts. We still assume all the 
eigenvalues to be distinct. Since L is a real matrix, the complex eigenvalues 
appear as pairs. For a convenience, we also assume that there is at most one 
pair having the same real part. Suppose A^ + i(3 k and A& — i(3 k are the first 
pair of complex eigenvalues. Then from ( |4.2j ), the leading order term in D k 
is 



#(A0 ... 0?(A fc + z/? fc ) 
0O(Ax) ... <f k (\ k + if3 k ) 



+ 
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0?(Ax) ... <f>i(\ k -i/3 k ) 2 



0° fc (Ai) 



Since is real by Proposition 1, one can write the above as 




[Acos(2p k t) + B sin(2/3 fc f)] . 



where A and B are two real constants. The above is an oscillating function 
about zero. Thus by Proposition 2, we conclude that L(t) blows up to infinity 
in finite time. □ 

This theorem implies that the complexness of the initial scattering data 
leads to the blowing up of the solutions. Now, only the situation left unde- 
termined is the case where S is indefinite, and the scattering data are real. 
We then have: 

Proposition 6 Suppose conditions as stated above are satisfied and detQ® ^ 
0, i = 1, • ■ ■ , N. If L(t) doesn't blow up to infinity in finite time, then L(t) 
has the sorting property, i.e., L(t) — > diag(X\, • • ■ , Xn) as t — > oo. 

Proof. If L(t) doesn't blow up to infinity in finite time, then Dj(t)'s are 
all positive definite, in particular, the leading order in the expansion Q4.2|) is 
positive. Then the same proof as in Theorem 3 applies here. □ 

5 Examples 

We here demonstrate our results using two simple examples. The first exam- 
ple includes a parameter, and we show a bifurcation behavior of the solutions 
as the parameter changes. The second example is for a blowing up aspect in 
the case where S is indefinite, and all eigenvalues of L are real and so is $°. 

Example 1. We take L to be a 2 x 2 matrix with S = diag(l, —1), that 
is, $ G 0(1, 1), and L is given by 




(5.1) 
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Our hierarchy (|2.1|) then gives for n=l 



dai _ _nu2 

dt ~ zo i> 
-26?, 



it 

da.2 
dt 



^ = -b 1 (a 1 + a 2 ) 



(5.2) 



From (|5.2| ), we can see both a\ and a 2 are always decreasing. If a\ + 02 is 
initially negative, then b\ increases faster and faster, and we expect a blowing 
up of the system. While, for the case where ai + a 2 is initially large positive, 
we expect that b\ to decrease to and a\ + a 2 decreases to some positive 
number. This is indeed the case. 

Let us take the initial conditions, a\ = 0, bi = 1, and a 2 = c, a constant. 
We then show a bifurcation behavior of the solutions with the parameter c, 
i.e., blowing up at c < c and the sorting at c > c for certain value of c . The 
initial scattering data, the eigenvalues Ai^ and the normalized eigenmatrix 
$°, are given by 



— c ± yc? 



(5.3) 



and 



a 2 



(5.4) 



The eigenvalues Ai and A2 take the following values as the function of c: > 
Ai > — 1 > A 2 , for c > 2; Ai = A 2 complex for |c| < 2; and Ai > 1 > A 2 > 
for c < —2. Note in particular that $° becomes complex (pure imaginary) 
when c < — 2, even though the eigenvalues are real. Then from (|3.17|) and 
( |3.21| ), we have 



D 1 (t) 



1 



A 2 — Ai 



(A 2 e 2Al * - A ie 



2X 2 t 



(5.5) 



and 



Aif 



_Ai__gA2* 



1 p A 2 t 

1 \A^! 



1 e Ait 



(5.6) 
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With the chioce of Ai and A2 in ( |5.3| ), we have Ai > A2 for all |c| > 2. In 
particular, notice that the dominant term in D(t) for large t becomes negative 
for the case c < — 2. This is due to the complexness of $°, and implies the 
blowing up in the solution at the time t = tg (Theorem 3), 

On the other hand, for c > 2 we have the sorting result. For |c| < 2, the 
eigenvalues become complex and Di(t) is expressed as 



D l {t) = -2(4 - c 2 )-*e- 2ct sin ( VA^t - 0) (5.8) 



with taxiO = V4 — ?/2. This also indicates the blowing up (Theorem 3). 



The solution L(t) is obtained from fl2.13|) 



~a n = ai =< Xtf >= A2g2Alf _ Aie2A2 , , (5.9) 

a 21 = h x =< A0 2 0! >= A2g2Alt _ Xie2X2t ■ (5-10) 

It is interesting to note that for the case c < — 2, the blowing up occurs 
at one time ts (|5.7p , and the solution L(t) will be sorted as t — > 00, with the 
asymptotic form, 

J _ x)- (5-H) 

Also note that if we start with a "wrong" ordering in the eigenvalues, i.e., 
Ai < A 2 , then we still have the correct sorting result with 

m (5.12) 

For the case of c = ±2, we obtain by taking the limit of ( |5.9|) and (|5.10| ) 

£(t)=f~# ^liV (5-13) 

\ l±2t \ 1* i±2t / 
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which showes the "sorting property" as t — > oo, i.e. L(t) — > ±diag(l, 1). It 
should be noted that L(0) at c = ±2 is not similar to ±diag(l, 1). 

In the above example, S is indefinite, when c > 2, the eigenvalues are real 
and so is $°, the system has the sorting property. The following example 
shows the blowing up aspect in that situation. 

Example 2. Take S = diag(l, —1,1), write 

L(0) = $°diag(2, 1, -100)S- 1 {$°) T S, 

where 




Calculating ( pT7|) , we get Di(t) = e 4t - 4e 2 * + 4e- 200 *. Since D 1 (0A) 



-3.95 < 0, we conclude L(t) blows up to infinity before t = 0.4. 



6 r-functions for the tridiagonal hierarchy 

In this section, we only consider the case of L being a tridiagonal matrix, 
and change t to 2t for convenience. L can be written as: 



/ siOi S2&1 

s\b\ s 2 a 2 s 3 b 2 

; 



V° 










s N-l a N-l S^bN-l 
s N~lbN-l SnO-N J 



(6.1) 



Similarly to the symmetric tridiagonal case, we can introduce r,'s as in ( |1.9|) . 
Equations (|1.9|), ( |1.1U| ) are now modified to, 



and 



s i s i+lt 



— log , 

Otl Ti-1 



r i+l T i-l 



n 



(6.2) 



(6.3) 
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However the equation for r» remains the same as ( 1.1 1|) . The r-functions are 
also given by the symmetric wroskian Ql.12 ) with g =< ((f) ( l(X)) 2 e^' t > >. 
Note in particular that the function g also satisfies ( |1.13[ ) with the change t 



to 2t, i.e. 



dg d n g 



(6.4) 



In the symmetric tridiagonal case, it was shown in 0] that r^'s are given by 
ZVs in (|3.17| ). In our case, Tj's are related to Dj's by 



Di 



Ti 



Si ■ ■ ■ Si 



(6.5) 



For instance, in example (i) of the previous section, T\ is given by D\ in 
(|5.5|) and T2 = —Di = — exp(— ct). It should be noted that in the symmetric 
tridiagonal case, these r-fuctions are positive definite, which is important in 
the moment problem of Hamburger, while in our case, the r-fuctions are no 
longer positive definite, in general. 

Let us now derive a hierarchy for Tj in Ql.12]) , which includes (|1.14j ) as its 
first member. First we note the following: In terms of P»(— <9 t ), where Pj are 
the elementary Schur polynomials defined by 



(6.6) 



n=l 



n-xn+i ■ 



Ql.llD can be written as 

where d t = (d tl , \d t2) ■ ■ •). As a generalization of ( |6.7D , we have: 
Lemma 4 For any k, I > 1 , the r-functions satisfy 



(6.7) 



Pk{-d t )Pi{-~d t )ri 



T; 



P k {-d t )n 



P fc -l(-4)P-l(-9t)T-,;_ 1 



Pii-d^n 



(6.8) 
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Proof. From the wronskian structure of the r$ with g satisfying 
have 



we 



Pk(-dt)T i = T i + 1 (i-k + l,i + l), (6.9) 

where 11+1(2 — k + l,i + 1) is the determinant of Tj+i after removing the 
i — k + 1th row and the i + 1th column. Note that Pk(— d t )ri = 0, for k > i. 
For 1 < k, I < z, we also get from the symmetric structure of the wronskian 

P k {-d t )Pi{-d t )n = r i+ i(i — k+l,i — l+ 1). (6.10) 

Noting Ti = r i+ i(i + 1, i + 1) we have 

P fc (-<9t)^(-4h] - [P k (-dt)%] [Pi(-dt)n 
= r i+l {i - k + l,i - I + 1)7*4.1 (i + l,i + l) 
-r i+1 (i — k + l,i + l)r i+1 (i + 1, i - I + 1). (6.11) 

The equation fl6.8|) then results directly from the the Jacobi formula for a 
determinant A{= r i+ i), 



A(i, k)A(j, I) - A(i, l)A(j, k) = A( l . k AA, 



(6.12) 



i k 



where A( ". ^ ) is the determinant of A after removing the zth and jth rows 
and the fcth and Zth columns. □ 

We then multiply \~ k fi~ l on the both sides of ( |6.8|) , and sum up k and / 
to get, from (|6.6|), 

exp(- Efc fcF^tJ exp(- Ei ]JA)t<] ^ 
exp(-E fe ^^JrJ [exp(-E/ ]Jr3c,)n 
exp(- Efc feir^ J exp(- Ei jfrdt^Tir-i] r l+1 . (6.13) 



Thus we obtain: 

Proposition 7 T/ie T-junctions satisfy a "bilinear identity" generating the 
relations of ( \6.%) , 



r .( t _ e [A] - e[fj])Ti(t) - Ti (t - e[X])n(t - e[//]) 
= ^r^-i(t - e[A] - e[^])r m (t), 



(6.14) 



where e[X) = (A 1 , ~A 2 , 
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We can now obtain a hierarchy written in the Hirota bilinear form ( |6.14j ). 
By setting \i = A, fl6.14j ) reads 



Ti (t - 2e[A])r i (t) - n(t - e[X))n(t - e[A]) 
= ^r,_ 1 (t-2e[A])r m (t) . 

Then changing t — e[A] to t, we can rewrite ( |6.15| ) as 

r f (t-e[A])r<(t + e[A]) -^(t)r,(t) 
= ^Ti^(t - e[X])r i+1 (t + e[X}) . 

This gives: 

Corollary 2 the Hirota derivatives D n := D tn in ( \L.ldQ , we have 

(exp(.D) - l) Tj • Ti = — exp(Z))r i+1 • r,_i , 



(6.15) 



(6.16) 



(6.17) 



where D:=E?7^D r 



Expanding (|6.17 ) in the power of A 1 , and using ( |6.6j ), we obtain equa- 
tions written in Hirota's bilinear form, for n > 0, 



P n+2 {D)Ti ■ n = P n (D)r i+1 ■ n-! 
The first three equations in ( |6.18| ) 



are 



(6.18) 



(6.19) 



(D 1 D 2 )r i ■ n = 2DxT i+ i ■ Tj_i 



(6.20) 



(Df + 3D 2 2 + 8D 1 D 3 )r i ■ n = 12{D 2 + Dl)r i+1 ■ . 



(6.21) 



Note here that ( |6.19| ) is just ( p.,14[ ) with t changes to 2t, and the "odd" 
powers of D^s on the l.h.s. of (|6.18|) are suppressed by operating on Tj • r, ; . 
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